Feedback control of the fluorescence light squeezing 
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We consider a two-level atom stimulated by a coherent monochromatic laser and we study how to 
enhance the squeezing of the fluorescence light and of the atom itself in the presence of a Wiseman- 
Milburn feedback mechanism, based on the homodyne detection of a fraction of the emitted light. 
Besides analyzing the effect of the control parameters on the squeezing properties of the light and of 
the atom, we also discuss the relations among these. The problem is tackled inside the framework 
of quantum trajectory theory. 
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Photo-detection theory in continuous time has been 
widely developed [H, 0, 0,13] and applied, in particular, to 
the fluorescence light emitted by a two-level atom stimu- 
lated by a coherent monochromatic laser [3, HI- As well as 
various feedback schemes on the atom evolution, based on 
the outcoming photocurrent, have been proposed @, 0|- 
However the introduction and the analysis of feedback 
have been mainly focused on the control of the atom Q • 

Here, we are interested not only in the atom, but also, 
and mainly, in the emitted light and in employing control 
and feedback processes to enhance the squeezing proper- 
ties of both systems. The squeezing of the fluorescence 
light can be checked by homodyne detection and spectral 
analysis of the output current [1, 0] . For these reasons we 
consider the mathematical description of photo-detection 
based on quantum trajectories, as it is suitable both to 
consistently compute the homodyne spectrum of fluores- 
cence light, and to introduce feedback and control in the 
mathematical formulation. We study how the squeezing 
depends on the various control parameters, how feedback 
mechanisms can be successfully introduced and which is 
the relationship among the squeezing properties of the 
quantum systems involved. We consider only Markovian 
feedback schemes a la Wiseman-Milburn [6| , as they leave 
the homodyne spectrum explicitly computable. 

Two-level atom and atomic squeezing. Consider a 
two-level atom with Hilbert space C 2 and lowering and 
rising operators cr_ and a + . Let us denote by a the vector 
{o~ x , o~y,o~ z ) of the Pauli matrices. Let the eigenprojectors 
of er z be denoted by P+ and P_ and, for every angle 0, 
let us introduce the unitary selfadjoint operator 

(70 = e 11 ^ <7_ + e~ 1( ^ cr + = cos 4><7 X + sin </> a y . 

A state p of the atom is represented by a point x in the 
Bloch sphere, p = (1 + x ■ a) /2, with fel 3 , \x\ < 1. 

Walls and Zoller [lOj suggested to define the squeezing 
of a two-level atom from the Heisenberg-Robertson un- 
certainty relations for a x and o~ y . By using the equatorial 
component a± of a with minimum variance Aer^ , we say 
that the atomic state p is squeezed if 



We call atomic squeezing parameter of p the quantity 

AS p = l-(x 2 +y 2 )-\z\, 

so that p is squeezed if AS P < 0. Thus min p AS P = — 1/4, 
attained by pure states with x 2 +y 2 = 3/4 and z 2 = 1/4. 

Detection and feedback scheme. We admit an open 
Markovian evolution for the atom, subjected to 'dephas- 
ing' effects and to interactions both with a thermal bath 
and with the electromagnetic field, via absorption and 
emission of photons. The atom is stimulated by a coher- 
ent monochromatic laser and the emitted light is partially 
lost in the forward channel and partially gathered in two 
side channels for homodyne detection. 

Let the free Hamiltonian of the atom be Ct>ocr 2 /2, with 
uiq > 0. Let the natural line-width of the atom be 7 > 0, 
let the intensities of the dephasing and thermal effects be 
fed > and n > 0, let the stimulating laser have frequency 
uj > 0. Finally, let the Rabi frequency be £1 > and let 
Alo = luo — ui denote the detuning. 

Let the fractions of light emitted in the forward and in 
the two side channels be |ao| 2 , |ai| 2 , 1 0:2 1 2 , respectively 
(|a | 2 + |ai| 2 + \a 2 \ 2 — 1); for k = 1,2, we can say 
that I CKfe | 2 is the efficiency of the detector k. Let the 
initial phase of the local oscillator in each detector be d}. , 
included in the parameter £ C by setting -d^ — argaj,. 
To change "dk means to change the measuring apparatus. 
Let the two homodyne photocurrents be I\ and I 2 . 

We introduce a feedback scheme a la Wiseman-Milburn 
based on I\ . Assuming instantaneous feedback, we mod- 
ify the amplitude of the stimulating laser by adding a 
term ge~ luJt Ii(t)/^/j proportional to It, with the same 
frequency w and with initial phase possibly different from 
that of the original laser. Let this phase difference be (p. 




Aa 2 ± = l-(x 2 +y 2 )<\{cr z )\ = \z\ 
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Then the atom has a Markovian evolution, whether we 
condition its state on continuous monitoring of the pho- 
tocurrents, or we do not. Let us call a priori state r\ t the 
unconditioned one and let us call a posteriori state p t 
the conditioned one. Of course rjt is the mean of pt- Let 
us write the evolution equations in the rotating frame, 
where they result to be time-homogeneous. Let us in- 
troduce first the parameters c = \g\ \ao\/y/j > 0, and 
Aui c = Auj + c 7 | a;i | cos($i — tp) € R. The a priori state 
rjt is governed by the master equation drj t = Crjt dt, where 
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+ 7 (n+l-M) \ *-P<t+- ^{P+ , p} 
+ y(ai a- - ica v ) p{a.\ er+ + ica v ) 
- |{(M 2 -2c|a 1 |sin(tf 1 -^))p + + c 2 , p] . 

The a posteriori state pt is governed by the non-linear 
stochastic master equation 

Apt = Cpt dt + y/jV[ai a- — ica v ]pt dWi(t) 

+ ^V[u2 0--]pt&W 2 {t), (1) 

where T>[a]p = ap+pa*— pTr [(a + a*)p] for every matrix 
a, and where W\ and W2 are two independent standard 
Wiener processes. The two homodyne photocurrents are 
given by the generalized stochastic processes 



J fe (t) = ^j\a k \ Tr [a# k p t ] + W k {t). 



(2) 



We suppose that |ao| is assigned by experimental con- 
straints and that the control parameters are O, Aw, $1, 
i?2, c, <p and, eventually, |ai| and \a 2 \- Of course, if c = 
there is no feedback action on the atom, so that its a pri- 
ori dynamics is independent of the measurement process, 
that is of «i, «2 and if. On the contrary, if c > 0, then 
the a priori dynamics is modified by the feedback loop 
and it depends also on a±, tp and c. 

In the Bloch sphere language L is an affine map. Let 
its linear part be given by the matrix ~A, where 
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Auj c — y(c\ai \ cos(i?i + tp) + c 2 sin 2^ , 
— Auj c — j(c\ai\ cos(i9i + tp) + c? sin 2tp\ 



7^- + ri + 2kd — 2c|ai| sin^i cosy? + 2c cos 95), 
7^1 + 2n- 2c|ai| sin(i?i - <p) + 2c 2 ). 



Apart from the exceptional case det A = 0, which occurs 
if and only if kd = n = 0, |ai| = 1, 2csin(i?i — y») = 1, 
fisini?! = 0, Auj = — 7CCOs("i?i — tp), the a priori 
dynamics has a unique stable stationary state rj cq = 



(1 



a) /2, which is asymptotically reached by for 



every initial preparation of the atom: 



-7(1 - 2c|ai| sin(i?i - tpfj A" 



Homodyne incoherent spectrum and fluorescence light 
squeezing. In a description of photo-detection based on 
quantum trajectories fl|, ©j the electromagnetic field 
has been traced out. Nevertheless, this description is 
fully consistent with a model which includes also a quan- 
tum description of the electromagnetic field and of its 
interaction with the atom and where Ii and I2 are the 
outputs of measurements performed just on the emitted 
light [jj. Therefore an analysis of the homodyne pho- 
tocurrents can reveal properties of the light detected in 
the corresponding channels. In order to study the squeez- 
ing properties of the light in channel k, the fundamental 
tool is the incoherent spectrum of I k [12| 
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It is the limit of the normalized variance of the Fourier 
transform of the photocurrent I k ; as usual, we call inco- 
herent the part of the spectrum due to the fluctuations 
of the output. The asymptotic behaviour of the atomic 
a priori state rj t ensures that the limit exists. It is a pos- 
itive even function of its real argument p which can be 
computed from equations (TTJ) and (J2J> by Ito calculus and 
by the full theory of quantum continual measurements, 
which can provide the first and second moments of Ii 
0,13]. Thus, for every initial state of the atom, we get 



Sk(p) = 1 + 2j\a k \ 2 s k 



A 



A 2 



(3) 



where s k and t k are the vectors in R 3 defined as 
s k = (cosz9 fc , sintffc, 0) , 
t k = Tr 



e 1 ** o~ _ ?7 cq + e r/eq 0+ - Tr[cr^ fc ?7 oq ] rj cq 
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Each spectrum Sk depends on p, fed, n, Q, Auj, a k , c and 
tp. Moreover, S 2 depends on at, too. In the case c = 
(no feedback) each dependence on tp. disappears and S2 
becomes independent of at; then, Si — S2 if cti — ot 2 - 
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When c = 0, for every /j, and the value of 5fc(/z; i?fc) 
is the variance of a quadrature of the light in channel k, 
the value of Ski^'&k + Tt/2) is the variance of the con- 
jugate quadrature and Heisenberg-type relations imply 
that 5fc(/i;^fc)5 fc (M;i?fc + tt/2) > 1 [12]. The light in 
channel k is in a squeezed state if the variance of one 
quadrature is below the standard quantum limit, that is 
if Sk{n) < 1 for some \x and 

When c > 0, only the light in channel 2 is potentially 
available for homodyne detection with arbitrary $2, as 
well as it could be employed for different uses. On the 
contrary, the light in channel 1 is a part of the feedback 
loop, it has to be detected, and a change of di implies a 
change of the atomic dynamics and, so, of the state of the 
emitted light itself, not only a change of the quadrature 
under consideration. Of course, the spectrum S% can be 
considered also in this case, but when Si(/i) < 1 one 
can speak only of 'in-loop squeezing'. Its meaning and 
possible usefulness are discussed by Wiseman 
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For each channel we can give a measure of the 'mean 
squeezing' of the light by introducing the quantity 



^ = 2^ 



[srWfc)-i] = ki 2 



When < for some the light in channel k is 

surely squeezed, but the spectrum can go below 1 even if 
nfc(z9fe) is positive. Moreover, we introduce the squeezing 
parameter of the state of the light in channel 2 



] 2 =infn 2 (i? 2 ) 
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(4) 



Control of squeezing. We are interested in the squeez- 
ing properties of the fluorescence light and of the atom. 
Regarding the fluorescence light, we can consider the 
squeezing of the light in the channels 1 and 2. Regarding 
the atom, we can consider the squeezing of the a posteri- 
ori state pt and of the a priori state rjt, and in particular 
of its limit ^ cq . Let us stress that the definition of atomic 
squeezing does not depend on the fact that we are work- 
ing in the rotating frame. Let us start by investigating 
the effect of the control parameters. 

Independently of the presence of the feedback loop, 
every time a parameter \atk\ vanishes, the corresponding 
photocurrent Ik reduces to a pure white noise (shot noise 
due to the local oscillator) with spectrum Sk = 1. 

The case c = 0. In this case the dependence of each 
spectrum Sk on the corresponding |afe| reduces to the 
explicit multiplication coefficient in Q. Therefore, when 
the control parameters Q and Aw give squeezed light in 
channel k, the lowering of Sk under the shot noise level is 
anyhow directly proportional to the fraction of emitted 
light gathered in that channel. 

For 17 = and n — there is no fluorescence light in 
the long run, so that each photocurrent Ik asymptotically 
reduces to a pure white noise and Sk = 1- 



FIG. 1: Channel 1 




FIG. 2: Channel 2 



For O = and n > there is no dependence on $k 
and Sk > 1. In this case there is only thermal light 
with carrier frequency o>o, while the local oscillator is at 
frequency to. The result are two temperature dependent 
Lorentzian peaks at [i = ±Aw. 

When fl > 0, Sk becomes ^-dependent and it can go 
below the shot noise level. This fact means that some 
negative correlation between the terms of the photocur- 
rent ([2]) has been developed. Some examples are plotted 
for both channels. All the figures are always in the case 
7 = 1, kd = n = and |ai| 2 = \a 2 \ 2 = 0.45, with the 
control parameters used first to fix the position of the 
minima of Sk and, then, to have the lowest minima. Fig- 
ures [T] and [2] show S\ and S2 (which are equal in this case) 
for minima in \l = (line 1, with Auj = 0, f2 = 0.2976, 
i9 k = -7r/2) and in fj, = ±2.5 (line 2, with Auj = 1.8195, 
Q = 1.7988, z? fc = -0.1438). 

One could also compare the homodyne spectrum with 
and without kd and n, thus verifying that the squeezing 
is very sensitive to any small perturbation. 

Regarding the atom, AS^ C depends only on kd, n, Q, 
Aoj. In the case fed = n = 0, the condition AS,, CQ < 



becomes < 20 2 < AAui 2 + -f 2 and the minimum value 
of AS,, oq is -1/8, reached for 6Q 2 = 4Au 2 + 7 2 . 
The case c > 0. 



The optimal squeezing in channel 
1 is always found for = and the feedback loop is 
very helpful, giving good visible minima of Si also when 
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|ai| is not close to 1. For example, FIG. Q] shows Si 
for minima in /i = (line 3, with Auj = 0, = 0, 
c = 0.2936, ip - <&! = tt/2) and in li = ±2.5 (line 4, 
with Auj = 2.5499, = 0, c = 0.3772, a = -1.3354, 
<^ = —0.0646). The utility of the feedback scheme can 
be appreciated by comparing line 1 with line 3 and line 
2 with line 4. 

If we are interested in the light emitted in channel 2 
and if |ai| and |at2 1 are assigned by some constraints, 
then the squeezing in channel 2 can be enhanced by a 
feedback scheme based on the photocurrent coming from 
channel 1, but the feedback performance is not as good 
as it can be for the squeezing in channel 1 itself. FIG. [2] 
shows S*2 for minima in p = (line 3, with Auj = 0, 
Q = 0.2698, i?i = tt/2, c = 0.0896, tp = 0, 2 = -tt/2) 
and in fx = ±2.5 (line 4, with Auj = 1.6920, Q = 1.9276, 
0i = 2.8168, c = 0.1326, tp = 1.2460, 2 = -0.0851). 

Anyway, if the only constraint is |ai| 2 +|a2| 2 = 1— |«o| 2 
and we are free in the choice of |ai| and |aa|, then the 
best observable squeezing in channel 2 is obtained in the 
case |ai| = 0, c = 0. That is, when the whole non-lost 
light is gathered just in channel 2 and the white noise I\ 
revealed in channel 1 is ignored. 

The feedback loop can be really efficient also to en- 
hance the atomic squeezing. For example, in the ideal 
situation |ai| = 1, fed = n = 0, with 7 = 1, Auj = 3, 
n = 4, 0x = tt/2, c = 1.3372, tp = -tt/40, we get 
AS^ = —0.2414, which is very close to the bound —1/4. 
In this case r] eq is almost pure so that also the a posteriori 
state p t is frozen in a neighbourhood of r] eq and AS pt is 
maximized, too. 

Fluorescence light vs atomic squeezing. There are not 
simple relations among the squeezing properties of fluo- 
rescence light in channel 1, of fluorescence light in chan- 
nel 2, of atomic a priori equilibrium state and of atomic a 
posteriori state. Indeed, changing the parameters of our 
model, we can observe a wide variety of behaviours. 

The only clear link is the one mentioned above: if 
AS^ ~ —1/4, then rj eq is almost pure and pt is frozen 
in a neighbourhood of ry cq , so that AS Pf is minimized, 
too, and the fluorescence light squeezing disappears as 
there is not incoherent scattering of light. One can check 
that actually only the coherent scattering survives, giv- 
ing a ^-contribution in p = to the complete spectrum. 
If the freezing of the atom is only approximate, one can 
check that all the spectra tend to become flatter and the 
fluorescence light squeezing tends to disappear. 

There is also the link (T5|) between AS^ oq and £2. This 
gives a direct relation between atomic and fluorescence 
light squeezing in absence of feedback. Indeed, in this 
case we have z cq < 0, so that £2 = 1^2 1 2 AS^ , and we 
can consider the case |ai| = 0, so that |ct2 1 2 is the fraction 
of the whole detected light. This relation is essentially 
the same found by Walls and Zoller considering a single 
mode for the emitted light [l(| ■ However the relation is 



not fundamental, as the feedback loop can give z cq > 
and in this case we have always £2 > even if AS^ cq < 0. 

There is no relation between fluorescence light squeez- 
ing revealed in channel 1 and in channel 2, even if we 
fix the constraint |ai| = |a 2 |. For example, the lowest 
minima of Si are found for c > and 17 = 0, but, every 
time ft = 0, the light in channel 2 is not squeezed as it 
can be proved that S2 > 1 for every p and every 2 . 

It is worth mentioning also the case \ot\\ — 1, 7 = 1, 
fc d = n = 0, with Auj = 0, Q = 0, 0i = tt/2, c = 1.2818, 
tp = 0. Then we have an extremely visible squeezing in 
channel 1 (Si reaches 0.3183), there is no squeezing of 
the atomic a priori equilibrium state (AS^ cq = 0.0922), 
while numerical simulations show that the a posteriori 
state pt tends to become pure (as |eti| = 1) with AS Pt 
stochastically moving between —1/4 and 0. 

Finally let us remark that the idea of the papers [7] is to 
choose the control parameters in such a way that, in the 
rotating frame, the atom is frozen in a preassigned pure 
state ho € C 2 , i.e. in such a way that, in the rotating 
frame, both the a priori state r\t and the a posteriori 
state p t asymptotically reach r\ cq — \h )(h \. This is 
possible in an exact way only in a very ideal case, which 
in our notations corresponds to |ai| = 1, fed = n = 0, 
Aw = 0, 0i = ±7r/2, tp = Q, which implies in particular 
012 = «2i = and x cq = 0. Then, p t is driven to a pure 
given state if f2 and c are such that y e 2 + z c 2 = 1 and 
2c sin 0i = 1 + z cq . But this implies t\ — and the two 
incoherent spectra reduce to pure shot noise. 
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